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^ . ABSTRACT 

The colour-magnitude relations of one of the most massive (w IO^^Mq), high red- 
\ shift [z — 0.89) clusters of galaxies known have been studied. Photometry has 

^ ■ been measured in the V, i?, /, z, F606W, F8I4W, J and K bands to a depth of 

' K « 20.5 = K* + 2.5 and spectroscopy confirms 27 K band selected galaxies as 

Tij" , members of the cluster. The V — K colours are equivalent to a rest-frame colour of 

T-H ■ 27OOA — J, and provide a very sensitive measure of star-formation activity. Hubble 

I Space Telescope imaging with the Advanced Camera for Surveys has been used to 

, morphologically classify the galaxies. 

The cluster has a low early-type fraction compared to nearby clusters, with only 
33 per cent of the confirmed cluster members having types E or SO. 
^ ■ The early-type member galaxies form a clear red-sequence in all colours. The 

Qh' scatter and slope of the relations show no evolution compared to the equivalent Coma 

Q , cluster relations, suggesting the stellar populations are already very old at z = 0.89. 

^ ■ The normalisation of the relations has been compared to models based on synthetic 

^ I stellar populations, and are most consistent with stellar populations forming at Zf ^ 3. 

. Some galaxies of late-type morphology were found to lie on the red-sequence of the 

ILJ ' colour-magnitude relation, suggesting that they have very similar stellar populations 

to the early-type galaxies. 

These results present a picture of a cluster in which the early-type galaxies are 
' all old, but in which there must be future morphological transformation of galaxies to 

match the early-type fraction of nearby clusters. In order to preserve the tight colour- 
magnitude relation of early-types seen in nearby clusters, the late-type galaxies must 
transform their colours, through the cessation of star-formation, before the morpho- 
logical transformation occurs. Such evolution is observed in the late-types lying on the 
colour-magnitude relation. 
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1 INTRODUCTION 

Observations of galaxies in clusters reveal several different 
evolutionary trends. On the one hand there is an apparent 
change in the galaxy populations of clusters with redshift. 
Whilst the cores of nearby clusters are dominated by early- 
type elliptical (E) and lenticular (SO) galaxies jPresslerl 
Il980r) . at high redshift there exists a significant population of 
blue, actively star-forming and po st-starburst galaxies whic h 
are almost absent at low redshifts jButcher fc OemleillQS^^ ■ 
The higher blue fraction of galaxies at high redshift may in- 
dicate that the morphological composition of the clusters is 
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evolving, with high redshift clus ters containing a larger frac- 
tion of spiral galaxies. Indeed, iDressler et al.l l|^997) show 
that there is a decrease in the fraction of SOs with increas- 
ing redshift, and a corresponding increase in the fraction 
of spirals. These observations seem to require some trans- 
formation of spirals to SOs, coinciding with a decrease in 
star- formation rate (e.g. iLarson et aLlTToSOll . possibly aris- 
ing from the infall of field spirals into the cluster and the 
subsequent quenching of the star-formation. 

On the other hand observations of early-type galaxies in 
clusters out to redshifts z ^ 1 reveal consistently old stellar 
populations, with correspondingly early assembly epochs. 
In particular, the tight correlation between th e colour and 
magnitudes of cluster galaxies at both low (e.g. lBower et all 
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| 1998[) and high (e.g. IStanford etaP Il998l : iHolden et aD 
120041) redshifts, the evolutio n of the fundamental plane 
fe.g. Ivan Dokkum et alJll99 8^ and the evolution of the K 
band luminosity function (e.g. Trentha^ fc Mobashor 1 99^ 
De Propris et alJ Il999l iNakata e"t"a"irboOll: iNelson eralJ 
200l| :|Kodama fc Boweil2003HToft et al.l2003l:lEllis fc Jonej 
2004l:lToft et al.l2004l:IStrazzullo et alj2006l) all reveal prop- 
erties consistent with formation at early epochs {z ^ 2), 
although LPs derived at bluer wavelengths shows a small 
fract ion of the stellar mass may be formed at lower redshifts 
(e.g. iNelson et al]|200ll : lAndreon et"al]|2004ri . 

Note also, that there are some problems with the inter- 
pretation of the Butcher-Oeml er effect as an i ndicat ion of 
morphological transformation. lAndreon et alJ (^93) show 
that the luminosity function of SO galaxies in the z = 0.023 
Coma cluster is consistent with the SO LF in the z = 0.41 
clus ter, which is at odds w ith an evolving population of 
SOs. iDe Propris et all i2003l) find very little increase in the 
blue fraction of galaxies when selection is made in the 
near-infrared, suggesting that the blue, high redshift galax- 
ies may be low mass systems rather than massive spirals. 
iBurstein et alj (|200^ show that SOs are brighter, and by 
implication more massive, on average than spirals in the 
NIR, which is inconsistent with a formation resulting from 
the stripping of mass from spiral galaxies. 

A further complication to the above picture is that nu- 
merical simulations of the hierarchical formation of clusters 
of galaxies reveal that as many as 50% of present day clus- 
ter galaxies may not have been within the cluster at z = 1, 
having s ince fallen into the cl uster as part of the assembly 
process jPe Lucia et al.l|2004^ . Thus any infalling galaxies 
which are progenitors for today's early-type galaxies must 
evolve to have remarkably similar properties to the galax- 
ies already in the cluster core, to account for the extremely 
small scatter of the local CMR of early- type galaxies. 

These issues have been studied in the CMRs of one of 
the most massive (M = 1.4 ± 0.2 X 10^^ Mo) cluster of 
galaxies known at z > 0.6 {z — 0.89). It thus provides 
an excellent opportunity to study the properties of galax- 
ies in high-density environments at high redshifts. Galaxies 
were selected in the near-infrared, providing a robust mea- 
surement of the luminosities of the old stellar populations, 
which is closely related to the dynamical mass of the galaxies 
((Gayazzi ct al. 1996). Spectroscopy provides accurate deter- 
mination of cluster membership for the brighter galaxies, 
and Hubble Space Telescope (HST) imaging with the Ad- 
vanced Camera for Surveys (ACS) provides morphological 
classifications. We have obtained near-infrared and optical 
photometry {V, R, I, z, F606W, F8I4W, J and K) from 
which CMRs have been constructed. The evolution of the 
galaxy population has been investigated via their scatter, 
normalisation and slope. 

The paper is organised as follows. In section (21 we de- 
scribe the cluster, and the sample. Section |2l describes the 
photometry. Section |1] details the morphological classifica- 
tion of the galaxies. In Section Owe discuss the fraction of 
early-type galaxies in the cluster. Section |^ describes the 
calculation of the CMR and presents the results. Section [7| 
investigates the evolution of the CMR via its scatter, nor- 
malisation and slope. Finally a discussion of the results and 
their implications for galaxy evolution in clusters are pre- 




Figure 1. UKIRT UFTI K band image of the z = 0.89 cluster 
C1J1226 with adapti vely smoothed XMM X-ray contours from 
iMaughan et alJl2004 The image is ~ 2.9 arcmin square. 



sented in section |H] A cosmology of Ho ~ 70km s ^ Mpc 
SIm = 0.3 and f^A = 0.7 has been assumed throughout. 



2 CLJ1226. 9+3332 

C1J1226.9-H3332, hereafter C1J1226 jEbeling et al.l I2OOII : 
ICagnoni etHlbOOll : iMaughan et alJl2004^ is at a redshift 
z = 0.892 ± 0.007. It has a bolometric X-ray luminosity 
of Lx ~ 5.3lo'2 X 10*^ ergs s~^, an X-ray temperature of 
Tx = ll.Sln g keV and a to tal mass Mtotai = l-4lo:2 x 10^^ 
Mq iMaughan et alj |2004|) . making it one of the most 
massive clusters known at high redshifts. The cluster ap- 
pears remarkably relaxed in X-rays indicating little dynamic 
activity, although recent XMM observations reveal tem- 
perature substructure and entropy profiles which suggest 
that the cluster may be undergoing a line-of-sight merger 
(Ben Maughan, private communication). The K band im- 
age with adaptively smoothed XMM X-ray contours (from 
IMaughan et alJl2004^ is shown in Figure Q 



2.1 Observations 

We have obtained near-infrared imaging and photometry of 
C1J1226, from which the initial galaxy selection was made, 
and a study of the K band galaxy luminosity function was 
made ( EUis,^ Jones 2004). This was followed up with opti- 
cal imaging and photometry, including high spatial resolu- 
tion HST imaging from which morphological classifications 
were made, as well as multi-object spectroscopy. These ob- 
servations will now be described in turn. 
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2.1.1 Near- Infrared imaging and photometry 

Observations were made on the 18th, 19th and 20th April 
2001 at the 3.8m United Kingdom Infrared Telescope 
(UKIRT) on Mauna Kea, Hawaii using the UKIRT Fast 
Track Imager (UFTI, Roche 200i) camera, a 1024 x 1024 
pixel HgCdTe array, with a pixel size of 0.091 arcseconds. 
Observations were made with the K98 and J98 filters (50% 
cut-offs 2.03-2.37/im and 1.17-1.33/im respectively) with an 
exposure time of 66s per frame and the seeing was typically 
« 0.5 arcsec in K. The total exposure times were 23946s in K 
and 26738s in J. The photometry was complete to K ~ 20.5 
mag and J « 21 mag, determined by the magnitude at which 
the number counts of galaxies begins to fall. A spatial dither- 
ing of 15 arcsec was employed and a mosaic was constructed 
from four pointings, since the field of view of 90x90 arcsec 
covered only a fracti on of the virial radius . On each night 
standard stars (from iHawarden et ahllioOll) were observed 
in order to calibrate the zero-point, air mass coefficient, and 
colour coefficient. One night had patchy cloud cover, and 
thus yielded inaccurate photometric calibration. This night 
was deemed to be non-photometric, and a calibration was 
obtained via observations of the same fields on photometric 
nights. Non-photometric data were then weighted according 
to their shorter effective exposure times. The observations 
covered an area of 2'47" x 3'1", corresponding to a physical 
size of 1.3 X 1.4 Mpc. 

2.1.2 Optical imaging and photometry 

We have obtained V, R, I and z band imaging of C1J1226 
at the 8.2m Subaru Telescope on Mauna Kea, Hawaii, using 
SuprimeCam, a prime focus camera consisting of 10 MIT/LL 
2048x4096 pixel arrays, with 0.2 arcseconds per pixel, yield- 
ing a total field of view of 34x27 arcminutes. 

The total exposure times were 2160s in V, 2880s in R, 
1920s in / and 1080s in z, yielding completeness limits of 
24, 24 , 24, 23 mag respectively. Standard stars from lLandold 
l|l992h were observed on each night. 

2.1.3 Multi-object spectroscopy 

In order to assess cluster membership, spectra were obtained 
and redshifts measured using the 10m Keck Telescope on 
Mauna Kea, with the Low Resolution Imaging Spectrom- 
eter (LRIS), and later using the 8.1m Gemini North Tele- 
scope, also on Mauna Kea, with the Gemini Multi-object 
spectrograph (GMOS). 

Forty-five galaxies have been confirmed as cluster mem- 
bers, 27 of which are within the UFTI K band field of 
view. The spectroscopic completeness of the K band selected 
galaxies is 82 per cent at isT < 18 and 61 per cent a.t K < 19. 

2.1.4 HST imaging 

As part of a larger project to study the lensing and morpho- 
logical properties of the galaxies in C1J1226 observations 
have been obtained with the 2.4m HST using the ACS. The 
ACS was used with the Wide Field Camera (WFC), which 
consists of two 2048 x 4096 pixel Scientific Imaging Tech- 
nology CCDs, with a field of view of 202 x 202 arcseconds. 
Six tiles of 4000s exposures each were mosaiced, in both 



the F606W and F8I4W filters, giving completeness limits 
of 25 and 23 mag respectively. 



3 PHOTOMETRY. 
3.1 Standard stars. 

The photometry of the galaxies was calibrated using stan- 
dard stars observed on the same nights as the cluster fields. 
The true magnitude of any celestial body is taken to be 

mtruo ^ZP- 2.51og(Ci?) + Asecz + B(J -~ K)t,uo (1) 

where ZP is the zero-point magnitude, CR is the count 
rate, A is the coefficient of extinction per unit airmass, z 
is the zenith angle {secz therefore being the airmass) and 
B is the colour coefficient necessary due to differences in 
the response of the combination of camera and filters used 
and the standard star system. Galactic extinction (0.06 mag 
in V, 0.007 mag in K), was assumed to be negligible given 
the possible systematic uncertainties of ~ 0.1 mag for each 
band. 

3.1.1 NIR, UKIRT 

Values of ZP, A and B were determined from observations of 
standard stars taken on each night. Standa rd stars were se- 
lected from the UKIRT faint standards list jHawarden et alJ 
I2OOII) . There were typically 6 stars of varying colours ob- 
served each night between 1 and 3 times, each at differing 
airmasses. 

Whilst the value of A and ZP may vary from night to 
night, the value of the colour coefficient B should be almost 
constant. Thus for the NIR a single value of B was mea- 
sured from all the photometric nights combined. The value 
we obtained is consistent with zer o, as expected from the 
filter design jTokunaea et al.ll2002^ . 

An independent value of ZP was measured for each 
night, and for the two photometric nights, an independent 
value of A was measured, and checked for consistency. Where 
this was not possible, on the cloudy night, a value of A 
from other nights was used and checked for consistency. All 
were found to be consistent, and a single value from one the 
photometric nights was used subsequently. In any case the 
low extinction in the K band and the low airmasses at which 
most observations were performed combine to make this a 
relatively small correction (<0.2 mag). 

Note that the final photometric error, as estimated from 
the scatter in the standard star measurements, is less than 
0.1 mag. The photometric calibrations are summarised in 
Table El 

3.1.2 Optical, Subaru 

For the optical SuprimeCam data standard star fields from 
iLandohi (£992) were observed each night. There were typ- 
ically two fields per night per filter, viz. SAlOl, SA98 and 
SA95 for /, SA95 and SA92 for R and SA113 for V. Note 
that each one of these fields has at least 10 standard stars. 
The z band was calibrated from the Landolt V, B, R and / 
band photometry, using the observed c olour transformations 
given in table 7 of ISmith et"al] J2002ll . 
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It was not possible to constrain the airmass coefficients 
from the standard star observations so average extinction 
factors from Subaru were used. The colour coefficient, B, 
was poorly constrained in all bands, hence we have excluded 
this term in the calibration. 

The zero-point values for each night were consistent, 
with final photometric errors of < 0.05 mag in all bands. 
The photometric calibrations are summarised in tabled 

3.1.3 Optical, HST 

The optical ACS data was calibrated using the known cal- 
ibrations provided in the fits image headers, and available 
from the Space Telescope Science Institute web-site^. The 
magnitudes are computed such that Vega would have zero 
magnitude in all bands. 

3.2 Galaxies 

Object detection was made in the K band image. The NIR 
is closely correlated with the dynamical mass of a galaxy 
iCavazzi et al and therefore a good choice for select- 

ing galaxies in an unbiased way for evolution studies. 

The SExTRACTOR software of Bortin & Arnoutj (Il996l') 
was used to search for objects in each field. To detect ob- 
jects a threshold value per pixel must be chosen along with 
a minimum number of connected pixels. Because the final 
K band images are constructed from a jittered pattern of 
images, the depth of observation varied across the final im- 
age, being at its deepest in the centres and shallowest at the 
edges. Therefore detection of objects was done in two re- 
gions for each image, using different detection parameters, 
in order to reach as deep as possible in the centre of the im- 
age whilst avoiding spurious detections at the edges of the 
image. The minimum significance of the detections was ~ 4(t 
in the centres of the images and « 7a" at the edges, where 
a is the background RMS determined from counts over the 
whole image (and is thus an overestimate of a in the central 
region and an underestimate in the outer region). Faint pix- 
els surrounding deblended objects are assigned to one of the 
sub-objects with a probability based on the expected contri- 
bution at that pixel from each of the deblended objects (see 
iBertin fc Arnout3ll996l) . Reliability of the object detection, 
and in particular the handling of overlapping objects, was 
checked by eye. 

Counts were measured using an adaptive aperture based 
on Kron's algorithm (JKron 1980), and also in a circular aper- 
ture with radius chosen to maximise the signal to noise. 
For each object a value of stellarity was also measured (see 
iBertin fc Arnout j|l996l) . 

To determine the magnitudes of the objects equation 
was used with the following complication. It is unknown to 
start with what the true colours of the objects are, therefore 
the last term in equation cannot be determined. To cir- 
cumvent this problem an approximate magnitude was mea- 
sured in each filter neglecting the colour term in equation 
allowing an approximate colour to be determined. The ap- 
proximate colour is then multiplied by a correction factor 

^ http:/ /www. stsci.edu/hst/acs/analysis/zeropoints 



previously determined using the same technique to mea- 
sure the approximate colours of standard stars and their 
true colours. The average colour correction for the NIR was 
(.i^rL = 0-982. 

Note that for the determination of colours, magnitudes 
were derived from fixed, circular apertures to avoid the ef- 
fects of internal colour gradients, whereas the adaptive aper- 
ture magnitudes were used to derive pseudo-total magni- 
tudes. The use of adaptive aperture magnitudes using SEx- 
tr actor's K ron radius should avoid the problem described 
bv lAndreonI |2002) of underestimating fiuxes for galaxies 
with low central surface brightness. The overall reliability 
of the K band photometry was checked by comparing the 
field galaxy counts d erived from the offse t fields with deeper 
published results (see lEllis fc Jonesl2004ll . There is generally 
good agreement down to our limiting magnitudes. 

Because the J band image has the same pixel scale as 
the K band image it was possible to extract the galaxies at 
the positions measured in the K band image using SEx- 
tractor's double image mode. The galaxy colours were 
measured in the same size apertures as for the K band data. 

The V, R, I, z, F606W and F8I4 W data covered a much 
larger area and had a different pixel scale and thus it was not 
possible to utilise the double image mode of SExtractor 
to extract the galaxies. Therefore catalogues of all the galax- 
ies were made separately for each of the V, R, I, z, F606W 
and F8I4.W band images using the same size apertures for 
the colours as for the NIR. K band selected galaxies were 
then extracted from this list using iraf's xyxymatch task. 
This takes the positions of three objects which are common 
to two images (preferably bright stars) and computes the 
transformation needed in terms of linear shifts, magnifica- 
tions and rotations. Thus a catalogue of K selected galaxies 
can be collated in each band. 

A few objects detected in the K band were not detected 
in the optical bands due to either the faintness of the optical 
luminosity or the failure to deblend two very close K band 
sources due to the different extents of the galaxy profiles 
in different bands. Note however, that the CMRs are de- 
termined only for galaxies with if < 19 (except for fainter 
confirmed cluster members) and there is only one object not 
in F606W, F8I4 W, R and / with K < 19 and one more ob- 
ject in V. Both of these are due to the very small separation 
from another object and the failure to deblend and subse- 
quently match these in the visible bands. Thus at K < 19 
the matching of galaxies is fairly robust and will not make 
a large difference to the derived CMRs. 

Star-galaxy discrimination was determined using SEx- 
tractor'S stellarity parameter, in the K band. A cut-off of 
0.8 was selected to delimit the two classes of objects with 
those objects with values greater than 0.8 being excluded as 
stars. This value was confirmed by examination of the radial 
profiles of detected objects. It was found that objects with 
a stellarity greater than 0.8 had an almost constant Gaus- 
sian FWHM close to the value of the seeing, whereas objects 
with stellarity less than 0.8 had more extended profiles. The 
results are insensitive to the precise value of this cut-off as 
there were very few stars in each field. 
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Table 1. Summary of photometric calibration. 

Waveband ZP A B Colour for B 



K 22.37 ±0.01 -0.18 ±0.05 0.00 ± 0.09 J - K 

J 22.78 ±0.01 -0.03 ±0.03 -0.03 ± 0.03 J-K 

z 27.08 ± 0.02 (fixed) 

/ 26.26 ± 0.04 -0.03 (fixed) 

R 27.50 ± 0.03 -0.09 (fixed) 

V 27.22 ±0.03 -0.12 (fixed) 



4 MORPHOLOGICAL CLASSIFICATION 

HST ACS data have been used to morphologically classify 
the confirmed member galaxies. These data have a point 
spread function (PSF) of < 0.1 arcseconds and a pixel scale 
of 0.05 arcseconds. Morphological classifications were car- 
ried out using 2-dimensional surface brightness fitting and 
isophotal classifications, as well as by inspection of the im- 
ages by eye. We attemped to decompose the galaxies into 
bulge and disc components by fitting a lSersid (119681) profile 
to the bulge and an exponen tial profile to the disc follow- 
ing |Khosroshahi_et_ai] teOOd) . A pure elliptical galaxy (no 
disc) is judged by a very small disc-to-total luminosity ratio 
(D/T < 0.05) while a pure disc galaxy by D jT > 0.95. The 
quality of the fit is checked by the value of the reduced 
and visual inspection of the residuals. We do not rely on 
the 2D surface brightness fit results if the reduced y^> 2. 
We use the HST PSF simulator, Tiny Tim jKrislil 19951) to 
simulate the PSF at the position of each galaxy in the ACS 
image, although the fitting results are found to be robust to 
the variations of the PSF across the camera. The images in 
F8I4W band are used for this analysis 

In addition to the above method we also performed an 
isophotoal analysis of the member galaxies via inspection of 
the radial surface brightness and ellipticity profiles in two 
different bands {F8I4W and F606W). The final categorisa- 
tion into early or late type galaxies has been decided by tak- 
ing into account all the information from each technique on 
an individual basis. All galaxies later than SO were classified 
as late-types. Generally galaxies which required a two com- 
ponent (bulge+disc) surface brightness profile, with the disc 
of equal or greater luminosity than the bulge were considered 
to be late types. Bulge dominated systems were considered 
to be early-type galaxies unless there was a significant spi- 
ral or twisting structure in the isophotes. Some galaxies (the 
BCG and two others) were not possible to classify into any 
conventional early or late type. 




Figure 2. A comparison of the early-type fraction in C1J1226, 
indicated by the large closed circle, with those pres ented by 
Ivan Dokkum et al ]i2OO0l) ( compr ising of ob servations bv |Dressleil 
Il98(l star: lAndreon et al.lll99'^ triangles; IPressler et al.lll997l_ 
open circles~lLubir ^^ami ^^ 
■2000l cross: .van Dokkum et 



closed sn nare: i FahHna.nt et a ^ 
I200(i open squares). 



galaxies at brighter limiting magnitudes increases the spec- 
troscopic completeness of the sample, and it is found that the 
early type fraction remains low in all cases. Table |21 shows 
the early-type fractions at a range of limiting magnitudes. 

The low early-type fraction is indicative of a 'mor- 
phological Butcher-Oemler effect', and agrees with previ- 



e.g. lAndreon et alJll997l; 


Dressier et al.lll997l: 


11993: Ivan Dokkum et al. 


|2000|). The earlv- 



type fract ion as a function of redshift was presented by 
van Dok kum et alJ ^2000) and we show that the early-type 
fraction of C1J1226 is in excellent agreement with this in 
Figure H 



5 EARLY TYPE FRACTION 

Of our twenty-seven confirmed cluster members the frac- 
tion of early-type galaxies is 33 per cent (9 early-types; 15 
late-types; 2 unclassified and the merging BCG). This low 
value could be partly due to a selection effect, since it is 
easier to obtain redshifts for spectra of late-types due to 
prominent emission lines. However, the galaxy selection for 
spectroscopy was made in the K band which is insensitive to 
the presence of young stellar populations within the galaxies, 
and would therefore provide an unbiassed selection criterion, 
and focussed on red galaxies roughly within the colour-range 
of early-type galaxies ai z — 0.89. Furthermore, selecting 



6 COLOUR-MAGNITUDE RELATIONS 

CMRs were fit for three different sets: all galaxies excluding 
confirmed non-members, confirmed members, and confirmed 
early-type members. Note that although the BCG shows 
classic signs of recent mergers, displaying multiple nuclei, it 
has been included in the early-type galaxies for the purpose 
of fitting the relations, since it clearly lies on the CMR in all 
cases and its extreme magnitude thus facilitates computing 
accurate slopes, normalisations and scatters for the CMRs. 

In order to fit the CMRs, the data were first divided 
into magnitude bins. Each bin then contains a distribution 
of colours. An estimate of the average value of this distri- 
bution was then estimated using the bi-weight estimator 
jBeers et al.lll990l) . where the average takes into account 
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Table 2. The morphological fractions of confirmed cluster members. 



Limiting 
magnitude 


Spectroscopic 
completeness/ per cent 


Percentage of 
early-types late-types 


others 


Total 
number 


K < 19.6 


45 


33 


56 


11 


27 


K <19 


61 


30 


57 


13 


23 


K < 18.5 


73 


29 


59 


12 


17 


X < 18 


82 


38 


50 


12 


8 



the skew of the distribution (due to a tight distribution of 
red galaxies associated with the cluster, and a population of 
bluer foreground galaxies). Once the average for each mag- 
nitude bin was determined a straight line was fit using the 
method of least squares. This straight line is the colour mag- 
nitude relation. The results are given in Table|21and are plot- 
ted in Figure |3| along with the best fitting relations for the 
confirmed early-type member galaxies. The diagonal dotted 
fine shows the approximate detection limit; galaxies redder 
than this will not have been detected. This was determined 
simply from the faintest reliable detection in each band. 

The early-type member galaxies form an obvious red se- 
quence in all cases. However, some non early-type galaxies 
are also clearly on the CMR. The BCG is a classic merging 
galaxy, showing double nuclei, and three very close compan- 
ions, however it clearly lies on the CMR. Another intriguing 
result is the presence of late-type galaxies on the CMR. We 
have conservatively adopted only E and SO galaxies as early- 
types, all spirals are considered to be late-types. The late- 
type galaxies on the CMR show a class of passive spirals, 
with the same photometric properties as elliptical galaxies, 
apparently consisting of similar populations of old stars with 
no obvious star-formation. 



7 EVOLUTION OF THE 

COLOUR-MAGNITUDE RELATIONS 

The evolution of the CMR has been investigated via the 
scatter, the slope and the normalisation. These investi- 
gations make extensive use of comparison to the nearby, 
z = 0.023, Coma cluster, for which the data were kindly 
provided by Roberto de Propris and are taken from Eisen- 
hardt et al. (2005, in preparation). 

When comparing the CMR a.t z — 0.89 with that at 
z = 0.023 the shift in observed waveband must be accounted 
for. Therefore in comparing the high z CMRs with those of 
Coma two different relations are used. The relations from 
the same filters are used to compare normalisations (sec- 
tion 17.211 as the models used take into account the effects 
of redshifts and fc-corrections. However, for comparison of 
the slopes and scatters the nearest equivalent Coma rela- 
tions are used, taking into account the shift in rest-frame 
wavebands when observing clusters at high redshifts. The 
equivalent wavebands used are listed in Table 2] Note that 
the a,t z = 0.89 the V band is equivalent to rest-frame wave- 
lengths of ~ 2700A. This is somewhat bluer than the U band 
used as the nearest equivalent Coma waveband, and thus we 
emphasise that caution must be used in deriving conclusions 
from the comparison of the V ~ K CMRs with the Coma 
U -J CMRs. 



7.1 Scatter 

The intrinsic scatter of the CMR is an extremely useful in- 
dicator of the star-formation and merger history of cluster 
galaxies. However, to accurately quantify the intrinsic scat- 
ter it is first necessary to properly account for the contri- 
bution to the overall scatter introduced by statistical un- 
certainties inherent in measuring g alaxy magnitud e s. We 
have followed the general method of IStanford et all lll998h 
to measure the intrinsic scatter and we give the details of 
the method below. We define the intrinsic scatter, 7s as. 



-^S = ^CTobs -(^phot. (2) 



where (Tobs is the overall scatter of the CMR as observed 
and (Tphot is the scatter due to errors in the photometry. 
The scatter and its error are d etermined using t he biweight 
location and scale estimators jBeers et alJll990l) of the dif- 
ferences between the best fitting relation and the data. 

It is paramount that aphot is determined accurately if 
the intrinsic scatter is to be meaningful. We have followed 
[Stanford ct al. ( 1998j) in our method of determining (Tphot 
via Monte-Carlo simulations. An artificial CMR was pro- 
duced in which the magnitude of each galaxy is held at the 
value measured for each galaxy. The colour of each galaxy 
was then set to a preliminary value determined from the best 
fitting relation. This colour was then allowed to vary by a 
random amount determined from a pseudo-random Gaus- 
sian distribution approporiate to the measured error on the 
colour. The scatter from this artificial relation was then cal- 
culated as before. This process was repeated 1000 times to 
accurately quantify the photometric scatter, aphot and its 
error. 

The intrinsic scatter was then calculated using equa- 
tion |5| This process was repeated on three subsets of the 
K band selected galaxies: all galaxies except confirmed non- 
members with K < 19, all confirmed members, and all con- 
firmed E and SO galaxies plus the BCG. The results are 
shown in Table 0] 

Comparing the early-type relations to the Coma rela- 
tions it can be seen that the high redshift relations are in 
good agreement with the nearest equi valent Coma relations . 
These results support those found bv lStanford et al] (Il998h 
for a larger number of clusters extending to similar redshifts, 
in which the scatter for all their clusters out to z ~ 0.9 is 
seen to b e consistent with the equivalent Coma relations. 
Similarly, lEUis et all il997ft find no increase in scatter out 
to 2: ~ 0.54. Note that the inclusion of late-type galaxies 
increases the intrinsic scatter measured. 
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Figure 3. Colour-magnitude relations for C1J1226. Filled circles indicate spectroscopically confirmed early-type cluster members, while 
the filled squares are the confirmed late- type members, filled triangles are cluster members with unclassified or unusual morphology, 
starred symbols indicate confirmed non-members, and open squares are of unknown membership. The solid line shows the best fit to the 
early type galaxies and the BCG as discussed in the text. The dashed line is the equivalent Coma relation normalised at K* = 18.05, 
see section [7.3l The dotted line shows the approximate detection limit in each colour. 
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Table 3. The best fitting CMRs as fitted by the bi-weight estimator of lBeers et al.l ll99(lh . The CMRs have been fit separately for all 
galaxies (excluding confirmed non-members) with K < 19 mag, confirmed members, and confirmed early- type members (plus the BCG). 

All Members Early-types 







Intercept 


Slope 


Intercept 


Slope 


Intercept 


Slope 


J-K 




2.4 ± 0.3 


-0.03 ±0.02 


3.1 ±0.8 


-0.07 ±0.04 


3.2 ±0.7 


-0.08 ±0.04 


z-K 




4.32 ± 0.09 


-0.057 ±0.005 


5.0 ±0.8 


-0.09 ±0.04 


4.9 ±0.8 


-0.09 ±0.04 


I - K 




5.4 ±0.2 


-0.10 ±0.01 


5.7 ±0.7 


-0.12 ±0.04 


5.6 ±0.8 


-0.11 ±0.04 


R~K 




7.5 ± 0.5 


-0.14 ±0.03 


7.3 ±0.8 


-0.13 ±0.04 


7.2 ±0.9 


-0.12 ±0.05 


V - K 




12 ± 1 


-0.32 ±0.07 


9±1 


-0.12 ±0.07 


9± 1 


-0.12 ±0.07 


F606W - 


K 


9.3 ± 1.0 


-0.21 ±0.06 


8.6 ±0.7 


-0.17 ±0.04 


9± 1 


-0.15 ±0.07 


F814W- 


K 


5.8 ± 0.6 


-0.14 ±0.04 


5.8 ±0.6 


-0.14 ±0.03 


5.9 ± 1.0 


-0.13 ±0.05 



Table 4. The intrinsic scatter of the CMR, for all galaxies with iC < 19 mag except confirmed non-members, all confirmed members 
and all confirmed early- type members. Also shown are the equivalent Coma relations for elliptical galaxies. 

Scatter Nearest equivalent Scatter of 







All 


Members 


Early-type members 


z = 


CMR 


Coma ellipticals 


J-K 




0.3 ± 0.4 


0.1 ± 0.1 


0.1 ±0.1 


R- 


J 


0.05 ± 0.04 


z-K 




0.3 ± 0.5 


0.1 ± 0.1 


0.05 ±0.07 


V - 


J 


0.06 ±0.05 


I -K 




0.4 ±0.6 


0.1 ± 0.1 


0.05 ±0.04 


B - 


J 


0.06 ± 0.04 


R-K 




0.6 ± 0.8 


0.2 ± 0.2 


0.08 ± 0.05 


U - 


J 


0.09 ± 0.07 


V - K 




1.0 ± 0.7 


0.2 ± 0.2 


0.13 ± 0.10 


U - 


J 


0.09 ± 0.07 


F606W 


- K 


0.9 ± 0.9 


0.2 ± 0.1 


0.2 ±0.2 


u- 


J 


0.09 ± 0.07 


F814W 


- K 


0.4 ±0.6 


0.1 ± 0.1 


0.09 ± 0.09 


B - 


J 


0.06 ± 0.04 



7.2 Normalisation 

We have used the normalisation of the CMRs to measure the 
epoch of star-formation assuming an instantaneous burst of 
star-formation followed by passive evolution. In theory it 
is possible to construct more complex models. In practice, 
however accurately modelling the star-formation histories of 
high redshift galaxies is well beyond the capabilities achiev- 
able through colour information of a single cluster alone. 
Therefore certain assumptions have been made to fix free 
parameters in the models. 



7.2.1 Constructing the models 

iBruzual fc Chariot! (l2003fl provide libraries of synthetic stel- 
lar populations (SSPs) and their spectral evolution at un- 
equal intervals of time from zero to twenty Gyr in the range 
3200A to 9500A with a resolution of SA. We have followed 
the recommendation of using the Padova 1994 evolution- 
ary track s . We assume the initial mass function (IMF) of 
IChabrieil i2003l) . Evolutionary tracks are provided with a 
range of metallicities denoted by Z, expressed as a fraction 
of the total mass, hence Zq — 0.02. Figure 0] shows an ex- 
ample of a passively evolving, solar metallicity spectral en- 
ergy distribution evolves from an age of IGyr to 9Gyr, the 
youngest spectra are the brightest. 

These SSPs assume an instantaneous burst of star- 
formation. More elaborate star-formation histories may be 
modelled by combining multiple SSPs in an appropriate 
manner. Here, however, we stick to the simple cases of purely 
passive evolution following an initial burst of star-formation 
at various redshifts, and a model in which no evolution oc- 
curs. 

The spectra must be transformed to appear as they 
would to an observer. A redshift of formation, Z{, was chosen, 
corresponding to f = in the extracted spectra. Age-redshift 




1000 1 0* 105 



X/ A 

Figure 4. Evolution of a passively evolving spectral energy distri- 
bution with Salpeter initial mass function and solar metallicity. 
The spectra range from an age of 1 Gyr (brightest) to 9 Gyr 
(dimmest) in steps of 1 Gyr. The spectra are displayed at a lower 
resolution than used in the modelling for clarity. 



relations were calculated assuming a cosmology with Ho=70 
km s"^ Mpc"\ = 0.3 and fl/^ = 0.7. Thus the redshift 
at which each extracted spectrum would be seen was calcu- 
lated. Once the spectra were redshifted to appear as they 
would to an observer, they were convolved with the trans- 
mission function of the filter through which the observations 
to which the models will be compared were made. The spec- 
tra were then integrated yielding the flux per solar mass of 
the model galaxy. The flux was then converted to magni- 
tudes, assuming the magnitude of Vega is zero in all bands. 
The flux of Vega was calculated by integrating a calibrated 
model spectral energy distribution under the appropriate 
filter transmission. 
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Thus the evolution of the model galaxy spectra, ex- 
pressed in terms of pseudo absolute magnitudes, was cal- 
culated at various wavelengths for redshifts of formation, 
zt = 1.5, 3 and 5. The magnitudes in different bands were 
then subtracted yielding the evolution in colour of the chosen 
model. All the steps described above were achieved through 
use of code provided with the SSP libraries. 



7.2.2 Results 

The models are shown in Figure In the observer's wave- 
length frame the colours become progressively bluer with 
age, since the galaxy is observed at smaller redshifts. The 
evolving models display bluer colours than the no evolution 
model, in which the stellar populations do not evolve but 
are merely redshifted and k-corrected appropriately. Fur- 
thermore, at fixed redshift younger stellar populations have 
bluer colours. 

The models are compared to the value of t he best fitting 
CMR evaluated at K* ( C1J1226, K* = 18.05. lEUis fc .ToneJ 
120041: Coma, JC*=10.9. IPe Prooris et ai]ll99St ) where the 
error on the C1J1226 point is given by. 



iV-2 



(3) 



where c and m, are the intercept and slope of the best 
fitting relation respectively, and A'^ is the number of points. 
The error on the Coma point is just the error on the in- 
tercept. The Coma colours for the HST ACS filters were 
calculated usi ng the ACS to Landolt system transforma- 
tions given bv lSirianni et all (1200511 for the E/SO template 
at 2 = 0.023. 

At the low redshift of Coma the models are indis- 
tinguishable from each other and the Coma value can be 
achieved by all the models. At 2; = 0.89, however, the mod- 
els diverge. With the exception of the z band all colours are 
consistent with formations at Zf « 3 and Zf ~ 5. The z band 
is not well fit by any relation, presumably due to systematic 
errors in the photometirc calibration. This is very likely a 
result of the transformations from t he V, B, R and / band 
magnitudes of the standard stars of iLandoIa (Il992l) to the 
z band, since the transformations used were derived from 
standard stars, and may not be applicable for the spectrum 
of galaxies at z = 0.89. Note, however, that though there 
may be a systematic offset in the zeropoint of the z band, 
the measured scatter and slope of the CMR will not be af- 
fected by this. 



7.3 Slope 

The slopes of the CMRs for C1J1226 have been compared 
with those of Coma. The slopes in each colour for C1J1226 
and those for the Coma relations and the nearest equivalent 
Coma relations (Eisenhardt et al. 2004, in preparation) are 
listed in Table |S] Once again we point out that the equiva- 
lence between the relations is only approximate and for the 
V — K relation there may be a significant difference between 
its true present day equivalent and the nearest equivalent, 
U — J, relation used. The evolution of the slopes is illus- 
trated in Figure 13 The figure shows the best fitting CMR 



compared to the equivalent Coma relations, which has been 
nor malised at K* = 18.0 5 . 

iKodama fc Arimotd ^1997^ and iKodama etlH lll998l) 
present two models of evolution of the CMR, one in which 
the relation is driven by variation in the metallicities of 
galaxies of different magnitude, and one in which the se- 
quence is driven by age. Their models are capable of pre- 
dicting the evolutio n of the slope of the CMR with red- 
shift. Figure 4 from IKodama et alj (^98) shows that the 
expected evolution for a metallicity driven sequence is a 
slight steepening of the slope towards higher redshift, in the 
observed waveband, due largely to the blue-shifting of the 
rest-frame band-passes. In the rest-frame the slope of the 
CMR will not change very much at all until ages younger 
than ~ 4Gyr, wh en it will become flatter du e to the changes 
in metallicities fKoda ma fc Arimotd Il997l . IKodama et alJ 
1998, Gladdcrs ct al. 199^. There is little difference between 
their passive evolution, Zf = 4.5 model and their no evolu- 
tion model. A comparison of our best fitting slopes with that 
figure would seem to indicate the our results are consistent 
with their passively evolving metallicity sequence model, but 
inconsistent with their age-sequence model for which the 
slope is expected to evolve very rapidly with redshift. Es- 
sentially the sa me trends are also produced in the hierar- 
chical models of iKauffmann fc Charloll (|l99?jl in which the 
largest elliptical galaxies are formed through the mergers of 
the largest, and therefore most metal rich, progenitor disc 
galaxies. 

The slopes of the C1J1226 CMRs are generally in agree- 
ment, within errors, with the equivalent rest-frame Coma 
relations. However, the C1J1226 relations are systematically 
slightly fiatter than the equivalent Coma relations. Fur- 
thermore, the difference in the slopes is greater for bluer 
colours. Note that here we are measuring something differ- 
ent from the expected ste epening towards higher redshifts in 
the IKodama et all il998l) models, as those models compare 
low redshift and high redshift CMRs in the same observed 
colours, not the same rest-frame colours. 



8 DISCUSSION 

Colour-magnitude relations of one of the most massive, high 
redshift clusters known have been determined. The early- 
type galaxies comprise only a small fraction of the cluster, 
with only 33 per cent of the confirmed member galaxies 
having type E or SO. This is indicative of a 'morphologi- 



previous results 


Andreon et alJll997 


1: iDressler et alJll997l: 


ICouch et al.lll998 


;lvan Dokkum et al. 


|2000|). The confirmed 



early- type members form a clear red-sequence in all colours. 
The early-type galaxies display remarkably similar proper- 
ties to early-type galaxies found in nearby clusters, suggest- 
ing that they are composed of similarly old stellar popula- 
tions. 

The scatter observed is consistent with the scatter in the 
nearest equivalent colours determined for elliptical galaxies 
in the Coma cluster within errors (Eisenhardt et al. 2004, in 
preparation), and suggests the stellar populations of early- 
type galaxies in C1J1226 are already old at 2; = 0.89. This 
is in agreement with previous wo rk which has fou nd no in- 
crease in scatter out to z ~ 0. 5 jEllis et aLlliggTi) . z « 0.9 




Figure 5. The normalisation of the OMR for C1J1226 (closed circle) and Coma (open circle) compared to the evolution of models of a 
simple stellar population with Z = Zq, Chabrier IMF, and Zf =1.5, 3, 5 respectively. Also shown is a no-evolution model. 
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Table 5. Comparison of the best fitting colour-magnitude relations for C1J1226 and Coma. 

C1J1226 Coma 

Equivalent observed-frame Equivalent rest-frame 



Colour 




Slope 


slope 


Colour 


slope 


J-K 




-0.08 ±0.04 


-0.024 ± 0.008 


R- 


J 


-0.088 ± 0.008 


z- K 




-0.09 ± 0.04 


-0.069 ±0.011 


V - 


J 


-0.101 ± 0.009 


I - K 




-0.11 ±0.04 


-0.092 ± 0.010 


B - 


J 


-0.146 ± 0.010 


R- K 




-0.12 ±0.05 


-0.113 ± 0.012 


U - 


J 


-0.220 ±0.014 


V - K 




-0.12 ±0.07 


-0.126 ± 0.013 


U - 


J 


-0.220 ±0.014 


F8I4W 


- K 


-0.13 ±0.05 


-0.092 ± 0.010 


B - 


J 


-0.146 ± 0.010 


F606W 


- K 


-0.15 ±0.07 


-0.113 ± 0.012 


U - 


J 


-0.220 ±0.014 



JStanford et alJ ll998D and even 2 ~ 1.3 iBlakeslee et alJ 
12001 iHolden et all2004) . lBower etHI lll99ai discuss the ef- 
fects of different star-formation histories on the scatter of the 
CMR. The basic principle is that star-formation spread over 
a larger time would result in a larger scatter. The colours 
also tend to approximately the same value > 5 Gyr after 
star- formation ceases. Thus the tight correlation of the CMR 
out to z = 0.89 is highly suggestive that the majority of the 
star-formation occurred within a relatively short space of 
time or muc h earlier than the ep och at which the galaxies 
are observed. iBower et alJ (|l99Sll also show that significant 
mergers would also disrupt the CMR causing an increase in 
scatter. Even if most gas is already in the form of stars, in 
which case very l ittle star-formation would b e triggered by 
a merging event jvan Dokkum fc Franxll20oil) . there would 
be an increase in the scatter of the CMR due to the mixing 
of galaxy colours. However, the effect is significantly smaller 
for for hierarchical merging as opposed to random merg- 
ing (iBower et al. 1998) . Kauffmann & Chariot (1998), show 
that the CMR can be reconstructed in a hierarchical sce- 
nario, if the progenitor galaxies have the necessary mass- 
metallicity relation. 

The slopes of the CMRs are also consistent with those of 
the nearest equival ent Coma relations. T his is in agreement 
with the results of IStan ford et al] (Il99^ and lHolden et al] 
i2004) . lK'odama et al.l (ll998^ show that the evolution of the 
slope of the CMR in a sample of 17 clusters is consistent with 
the relation being a metallicity sequence as opposed to an 
age sequence. Thus any processes which changes the metal- 
licity of a galaxy would disrupt the slope. A simple com- 
parison of the measured values with the evolutionary tracks 
of [Kodama ct al. ( 1998) suggest that our values are con- 
sistent with the Zi — 4.5, passively evolvin g metallicity se- 
quenc e presented. This is in agreement with lCladders et alJ 
il998l) who find an evolution of the slope of 6 clusters at 
0.2 < z < 0.75 and 44 clusters at z < 0.15 is consistent with 
models of lKodam a e^lj (^9^, al t hough direct comparison 
of our results with lCladdgr^t al.l (Il998l) is complicated by 
the different wavebands studied. 

Our result s are inconsistent with the findings of 
iFerreras fc Silkl i2000l) who find an increase in slope and 
scatter of the UV-optical CMR of CI 0939+4713 at z = 
0.41, which they argue is inconsistent with monolithic col- 
lapse. They also point out the sensitivity of using the 
UV part of the spectrum as an indicator of star for- 
mation. However, our results, do analyse the restframe 
UV part of the spectrum since at z = 0.89 V is sam- 
pling about 27OOA. The different results may be due to 
real differences in the star-formation histories of the dif- 



ferent clusters examined. C10939+4713 is thought to be 
part of a merging syste m iSchindler fc Wambsgansjll996l : 
iDe Fihppis et al.1 l2003^■ with a large fraction of post- 
starburst galaxies dBelloni et alj [l995). which may account 
for its enhanced star-formation. However, XMM-Newton ob- 
servations of C1J1226 (Ben Maughan, private communica- 
tion) also suggest that it is part of line-of-sight merger, but 
we do not see any evidence for enhanced star-formation. 
On the other hand the different r esults may be i n part 
due to differences in the analysis. IFerreras fc Silkl (l2000h 
use a sample base d on the morphological classifications of 
ISmail et alJ (|^93), which does not include redshift infor- 
mation, therefore the sample may be contaminated by fore- 
ground and ba c kgrou nd galaxies. The detection limits of 
IFerreras fc Silkl i2000l) may also affect the steep slope they 
measure, since they may exclude galaxies on the predicted 
CMR fainter than F702W 21.5, but fit to galaxies down 
to F702W ~ 23.5. Note, however, this would make the scat- 
ter they measure even larger, though again we point out that 
the CMR may include non-cluster members. Finally we note 
that we compu t e colo urs within a fixed aperture, whereas 
IFerreras fc Silkl (l2000h used the total magnitudes. 

This issue of aperture size is important in deriving the 
colour-magnitude relation. IScodeggia (1200 discusses the 
effect of using fixed aperture magnitudes compared to a fixed 
fraction of the galaxy light. They find that theU -V CMR 
of the Coma cluster is much flatter, and has a larger scatter 
when a flxed fraction of the galaxy light is used (a slope 
of -0.016 ±0.018 compared with -0.074 ± 0.008). They at- 
tribute these effects to the internal colour gradients of galax- 
ies, which have a large scatter, and are generally redder to- 
wards the centre of elliptical galaxies. Thus colours derived 
from fixed apertures will generally be redder for larger (and 
therefore usually brighter) ga laxies. However^ the v ery tight 
relations of the Coma CMR dTerlevich et alJl200lh derived 
from fixed size apertures are difficult to understand in this 
scheme, since the large scatter of the colour gradients would 
also produce a large scatter at fixed radius. An explanation 
may be that using a fixed aperture results in a higher sig- 
nal to noise, and hence a tighter relation, or perhaps if the 
colour gradients were generally small, and hence the scatter 
introduced would be small. Note that for the high-redshift 
galaxies reported here, the accurate computation of effec- 
tive radius, and accurate photometry in large apertures is 
not possible, and thus a fixed angular aperture was used. We 
investigated the effect of various aperture sizes on our results 
and found that for apertures larger than 1.5 arcseconds the 
measured colours were consistent. 

The effect of galaxy colour gradients could also ex- 
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plain the differe nce between the recent GALEX results from 
lYi et al.l (l2005l) . which show an increase in the scatter of 



the CMR when using near-ultraviolet— optical colours, de- 
rived from total magnitudes rather than aper ture magni- 
tudes. However, the res ults of IScodeggio i200j) cannot ex- 
plain the difference with lFerreras fc SiM ( 200(J) who also use 
total magnitudes, but measu re a steeper slope. 

The merging model of iKauffmann fc Charloli (Il998l) 
can reproduce the slope of the CMR at low redshifts due 
to the mass-metallicity relation of the progenitor galax- 
ies being preserved in the hierarchical nature of the merg- 
ing. These models also predict a progressive flattening of 
the CMR at high redshift, becoming significant at 2; > 
1 due to the younger average age of the most massive 
ell iptical galaxies. Some suppo rt for the merging model 
of IKauffmann fc Charlol lll998h is found in the BCG of 
C1J1226. High resolution ACS data reveals that the BCG 
displays classic signs of merging activity: double nuclei 
and several close companions. The BCG clearly lies on 
the CMR and has been included in measuring the evolu- 
tion. Such an obvious merger provides qualitative evidence 
that merging galaxies can remain on the CMR. However 
we caution that BCGs are not typical galaxies, and very 
likely have distinc t formation histories from the general 
cluste r pop ulation (iBhavsar fc Barrowlll985l) . lBrough et all 
J2002h and iBrough et al.l i200^ show that BCGs in clus- 
ters with high X-ray luminosity assembled most of their 
stellar mass at z > 1 and thus the subsequent assembly 
of the BCG is probably not amongst equal mass galax- 
ies, but is more likely a case of 'galactic cannibalism' in 
which the BCG is much more massive than the galaxies 
with which it is merging. Thus the merging activity may 
not have a la rge effect on the overa ll photometric properties 
of the BCG. lYamada et alJ (|20o3) report a similar merg- 
ing of the BCG in the z = 1 .26 cluster RXJ0848. 9-^4452, 
and lvan Dokkum et al.l (1200 jl show that the three brightest 
galaxies in the nearby cluster RXJ0848-I-4453 a,t z — 1.27 
have similarly undergone recent mergers. 

On the other hand the merging BCG may be dis- 
playing behaviour typical of many early-type galaxies in 
clusters. If mergers in clusters are predominantly between 
red galaxies, with old stell ar population, as in MS1054- 
03 jvan Dokkum et alJll999^ . the colour-magnitude relation 
may be only a little affected during the mergers, displaying 
an increased scatter, and wil l quickly come to resemble the 
CMR of passive early-types jTran et alJl2005ft . 

The normalisation of the CMR shows evolution in the 
sense that galaxies at high redshift are generally bluer than 
their present day counterparts (with the exception of z — K, 
for which there may be significant systematic uncertainties 
in the photometric calibration, and J — K which is consis- 
tent with no evolution within the errors), once the change in 
observed waveband has been taken into account. This evo- 
lution is consistent with purely passive evolution in all cases 
except for z — K, with a redshift of f ormation zf > 3 . These 
res ults are consistent w ith those of IStanford et al.l (Il998l) 
and lHolden et all (1200411 . who fit CMRs consistent with pas- 
sive evolution. 

Thus the elliptical galaxies seem to be composed pri- 
marily of old stellar populations with very similar properties 
to those found in nearby systems. This is consi stent with the 
evolution of the K band luminosity function jEUis fc Jonej 



l2004fl in which it is seen that the assembly history of the 
galaxies is consistent with monolithic collapse at zf ^ 2. 

An important point however, is that we have selected 
early-type galaxies to define our colour-magnitude rela- 
tion. Therefore we may be pre-selecting galaxies which con- 
tain only old stellar popul ations; the progenitor bias of 
Ivan Dokkum fc Franxl (|200lln . Note if this is the case, then it 
follows that any galaxies either outside clusters at z = 0.89, 
or of different morphological type, which are the progeni- 
tors of present day early-type galaxies in clusters, must end 
up with stellar populations indistinguishable from the early 
type galaxies already present in clusters at z — 0.89, since 
the CMRs of clusters at high and low redshift are so simi- 
lar. This could be the case if the progenitors are late-types 
containing similar populations of old-stellar populations, as 
well as younger populations. Once star-formation ceases, the 
galaxy will quickly become redder, as the massive, hot, lu- 
minous, short-lived blue stars turn off the main sequence 
and enter the giant phase. Thus they could rather quickly 
resemble the old stellar populations of early-types soon af- 
ter star-formation has ceased. Of course at some point there 
must also be some morphological transformation. 

The core C1J1226 seems to be showing such evolution 
in process. There are four-late-type galaxies within 0.1 mag- 
nitudes of the best fitting CMR inV — K, which is the most 
discriminating CMR due to its large scatter. These galax- 
ies all possess substantial discs. It would seem that these 
galaxies typify the process described above, whereby star- 
formation ceases in spirals and they migrate onto the CMR. 
Their presence on the CMR is suggestive that evolution of 
the stellar-populations occurs before any morphological evo- 
lution which may take place. This order of stellar evolution 
followed by any morphological evolution, is n ecessary to pre- 
se rve the CMR at all redshifts, as observed bv lStanford et all 
11998), otherwise the CMR would exhibit evolution in its 
scatter. 

The presence of passive spirals have previously been 
observed in high a nd intermediate redsh i ft clusters (e.g. 
ICouch et al.' '1998'; 'Dressier et alj Il999l iPoggianti et alJ 
il999c lYMnauchi fc Goto 2004 ) and in low redshi f t clus - 
ters ([van den Bergh 1 97(^ : llCoopmann fc KennevI Il99^ . 
iGoto et^lTl2003 ) observe that the passive spirals in these 
systems are found preferentially in the outskirts of clus- 
ters, where recently infalling galaxies are likely to be found, 
thereby providing a mechanism for the cessation of the star- 
formation, as the galaxies interact with the high-density en- 
vironment. Here however we find such passive galaxies in the 
core of an extremely massive, virialised cluster. The pres- 
ence of passive late-types in such a dense environment may 
be problematic if the transformation is entirely due to the 
infall of the galaxies. Thus it may be that age is also impor- 
tant in the cessation of star-formation. Galaxies in high den- 
sity environments are likely to be older, since the high over- 
densities would collapse at earlier epochs in an expanding 
universe. Therefore they could have already consumed their 
supply of gas by z = 0.89, and be in the process of trans- 
formation. Under such an explanation the morphology- and 
star-formation-density relation (eg. Lewis ct al. 2002) would 
be a refiection of an underlying age-density relation, due to 
the fact that galaxies form earlier on average in high-density 
environments. 

Stronger constraints on the degree of morphological evo- 
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lution and star-formation history should be possible via de- 
tailed spectroscopic and morphological studies of high red- 
shift clusters as a function of density and environment. The 
ACS data presented here will be further exploited in such 
studios to reveal finer details of the galaxy population of this 
remarkable cluster. 
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